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Preface 


The  purpose  of  this  study  was  to  measure  the  temporal- 
spatial  temperature  distributions  In  a  non-sootlng  CO-H2 
diffusion  flame  at  high  repetition  rates.  The  observations 
of  such  a  flame  at  Re  =  857  revealed  the  existence  of  both 
low  frequency  oscillations  associated  with  large  scale 
buoyancy-driven  vortices  and  relatively  high  frequency 
oscillations  believed  to  be  a  characteristic  of  the  Internal 
cold  let  fuel  core  of  the  flame.  Other  H2-CO:N2  gas 
mixtures  (using  the  H2  as  diluent)  were  Investigated  up  to 
Re  =  2264,  and  similar  high  frequency  oscillations  were 
found . 

For  this  research  the  novel  thin  filament  pyrometry 
(TFF)  technique,  developed  by  Goss,  Vllllmpoc  and  Sarka,  was 
employed;  this  technique  relies  upon  the  blackbody  emission 
of  a  thin  ceramic  filament  (SIC)  In  a  combusting  flowfleld. 
Radial  temperature  profiles  were  taken  by  the  TPF  technique 
In  the  H2-CO  flames  at  10  mm  axial  Increments. 

A  laser-sheet  lighting  technique  consisting  of  two 
perpendicular  sheets  was  used  to  visualise  the  flame 
structure:  this  technique  Involves  seeding  both  the  central 
fuel  Jet  and  co-annular  air  with  TiCl4  vapor.  _ 
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Abstract 

The  purpose  of  this  Investigation  was  to  measure  the 
radial  temperature  profiles  as  a  function  of  time  In  a  non- 
sootlng  CO-H2  diffusion  flame  at  high  repetition  rates. 
Measurements  were  taken  from  the  nozzle  exit  up  to  200  mm 
above  the  nozzle  In  10  mm  Increments.  The  results  of  such 
an  analysis  at  Re  =  857  revealed  the  existence  of  low 
frequency  (25-30  Hz)  oscillating  structures  In  the  outside 
reaction  layer  associated  with  the  dynamics  of  large  scale 
buoyancy-driven  vortices,  and  high  frequency  (300-325  Hz) 
oscillations  believed  to  be  associated  with  local 
Instabilities  of  the  Internal  Jet  core  of  the  flame. 

Other  H2-C0:N2  gas  mixtures  (at  Re  =  1717.  2099  and 
2264)  were  Investigated  and  similar  high  frequency 
oscillations  were  found. 

In  this  research  the  novel  thin  filament  pyrometry 
technique,  developed  by  Goss.  Vllllmpoc  and  Sarka  was 
employed:  this  technique  relies  upon  the  blackbody  emission 
of  a  thin  ceramic  filament  (SIC)  In  a  combusting  flowfleld. 

A  limited  Fourier  Transform  frequency  analysis  was 
carried  out.  However,  due  to  the  limited  number  of  sampled 
Intervals  an  exact  correlation  between  the  different 
harmonics  that  appeared  In  the  spectra  can  not  be  made. 
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EXPERIMENTAL  INVESTIGATION  OF  A  CO-H2  DIFFUSION 


FLAME  USING  HIGH-FREQUENCY  THERMOMETRY 


I. 


Background 

Combustion  processes  applied  to  propulsion,  energy 
conversion,  manufacturing  and  heating  are  vital  to  modern 
day  civilization.  Despite  the  century-long  use  of 
combustion,  there  are  a  number  of  areas  that  are  not 
currently  well  understood.  Many  of  the  current  limitations 
in  the  knowledge  base  of  combustion  can  be  attributed  to  the 
inability  to  experimentally  probe  these  processes.  These 
difficulties  arise  due  to  both  the  hostility  and 
interactivity  of  combustion  phenomena  which  makes  physical 
probing  quite  difficult.  Optical  techniques  have  long  been 
recognized  for  their  potential  to  non-intrusively 
Interrogate  combustion  processes.  With  the  advent  of 
lasers,  the  capabilities  of  optical  techniques  have  been 
greatly  enhanced.  Perhaps  more  importantly,  an  array  of  new 
methods  has  emerged  which  previously  was  not  possible  to 
Implement  or  was  yet  undiscovered.  Most  laser-based 
techniques  fall  in  this  latter  category  and  the  measurements 
of  combustion  temperature,  chemical  composition  and  flow 


velocity  are  frequently  based  upon  the  interpretation  of 
radiative  emission  spectra  which  follows  the  Interaction  of 
laser  light  with  the  specific  Internal  energy  states  of 
atoms  or  molecules. 

Since  combustion  problems  requiring  theoretical 
analysis  are  primarily  concerned  with  the  flow  of  reacting 
and  diffusing  gases,  one  must  understand  (in  addition  to 
chemical  thermodynamics)  the  conservation  equations  of  fluid 
dynamics,  including  transport  properties  and  chemical 
kinetics.  In  diffusion  flames.  Infinitely  fast  chemical 
reaction  rates  are  assumed  and  the  global  flame  properties 
are  determined  primarily  by  transport  effects. 

A  diffusion  flame  may  be  defined  as  any  flame  in  which 
the  fuel  and  oxidizer  initially  are  separated;  with  this 
usage,  the  term  is  synonymous  with  non-premixed  combustion. 
In  a  restricted  sense,  a  diffusion  flame  may  be  defined  as  a 
non-premixed.  quasi-steady,  nearly  Isobaric  flame  in  which 
most  of  the  reaction  occurs  in  a  narrow  zone  that  can  be 
approximated  as  a  surface  (flame  sheet  or  reaction-sheet). 

In  the  laminar  flame  regime,  all  the  mixing  and  transport  is 
achieved  by  molecular  processes,  while  in  the  turbulent 
flame  regime  these  processes  are  aided  by  macroscopic  eddy 
motion. 

The  most  fundamental  question  in  turbulent  combustion 
concerns  the  applicability  of  a  continuum  description.  It 
is  customary  to  assume  that  the  full  Navier-Stokes  equations 
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(that  is.  the  Navler-Stokes  equationa  augmented  by  chemloai- 
klnetlcs  laws)  govern  turbulent  combustion:  Navler-Stokes 
turbulence,  characterized  by  fluctuations  of  velocities  and 
state  variables  In  space  and  time,  occurs  when  the  Inertial 
forces  are  sufficiently  large  In  comparison  with  the  viscous 
forces  (that  Is.  at  sufficiently  large  values  of  a 
representative  Reynolds  number.  Re)  (1:374).  Turbulence 
develops  from  f luld-mechanlcal  local  Instabilities  of  flow 
(2:443-473)  when  Re  Is  large  enough  to  prevent  viscous 
damping  from  suppressing  fluctuations.  The  Navler-Stokes 
equations,  combined  with  a  proper  set  of  Initial  and 
boundary  conditions,  are  deterministic  In  that  they  are 
believed  to  possess  a  unique  solution.  However,  at  high  Re, 
the  solution  Is  strongly  sensitive  throughout  most  of  the 
flow  field  to  the  Initial  and  boundary  conditions,  and  In 
reality  the  conditions  cannot  be  specified  accurately  enough 
to  obtain  the  deterministic  solution,  either  theoretically 
or  experimentally.  It  then  becomes  reasonable  to  Introduce 
statistical  methods  and  seek  only  probabilistic  aspects  of 
the  solutions. 

For  Initially  non-premlxed  reactants,  two  limiting 
cases  may  be  visualized,  namely,  the  limit  In  which  the 
chemistry  Is  rapid  compared  with  the  fluid  mechanics  and  the 
limit  In  which  It  Is  slow.  The  primary  combustion  reactions 
In  ordinary  turbulent  diffusion  flames  encountered  In  the 
laboratory  and  In  the  Industry  appear  to  lie  closer  to  the 
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fast-chemistry  limit.  Some  of  the  specific  questions  asked 
are:  What  is  the  flame  length  or  flame  height?  What  is  the 
average  local  volumetric  rate  of  heat  release?  What  are  the 
initial  conditions  for  a  turbulent- Jet  diffusion  flame  to  be 
lifted  off  the  fuel  duct?  What  are  the  radiant  energy  fluxes 
emitted  by  the  flame?  How  complete  is  the  combustion?  What 
is  the  rate  of  production  of  oxides  of  nitrogen?  What  is  the 
rate  of  production  and  liberation  of  smoke  and  of  unburnt 
hydrocarbons? . 

In  order  to  test  the  theoretical  models  of  combustion, 
one  needs  to  be  able  to  measure  local  temperatures,  flow 
velocities,  and  concentrations  of  major  and  minor  species  at 
rates  that  are  in  excess  of  the  time  scales  of  turbulence. 

An  "ideal  "  technique  should  satisfy  the  following  criteria: 
It  should  (i)  be  capable  of  making  in  situ  and  real-time 
measurements,  (ii)  be  non-perturbing,  (ill)  have  a  high 
sensitivity,  (iv)  have  a  high  degree  of  spatial  resolution, 
and  (V)  have  a  high  degree  of  temporal  resolution.  In 
addition,  a  truly  ideal  technique  should  (vi)  measure 
several  combustion  parameters  simultaneously  and  (vii)  be 
affordable  (3:1339).  In  general,  optical  techniques  satisfy 
the  first  two  criteria  and  many  optical  techniques  have  been 
investigated  which  satisfy  the  criteria  enumerated  above  to 
varying  degrees.  A  few  of  the  techniques  have  been  quite 
successful . 

Temperature  and  majority  species  are  usually  determined 


by  means  of  Raman  techniques  (4:1521:  5:175).  while  radical 
Intermediates  In  a  much  lower  concentration  are  often  probed 
using  laser-induced  fluorescence  (LIF)  (6:568;  7:145: 

8:511).  Raman  methods,  however,  [with  the  exception  of 
electronically  enhanced  Raman  techniques  (9:121;  10:39)]  are 
Insensitive  to  trace  minority  species  since  their  signal 
strength  Is  proportional  to  the  concentration;  and  LIF 
techniques  often  suffer  from  colllslonal  quenching  of  the 
fluorescence  In  high-pressure  combustion  environments 
(8:511)  whose  effect  can  be  overcome  partially  by  using 
saturated-LIF  (11:445;  6:568). 

Newly  developed  techniques  Include  photoacoustic 
spectroscopy  (PAS)  (12:118;  13:2133).  photoacoustic 
deflection  spectroscopy  (PADS)  (14;  15:781),  and 
photothermal  deflection  spectroscopy  (PTDS)  (16:1333). 

After  Patel  and  Tam  (17:467)  demonstrated  In  1979  a 
highly  sensitive  pulsed  photoacoustic  spectroscopy  technique 
Involving  the  use  of  pulsed  lasers,  has  this  technique 
received  attention  In  the  area  of  combustion  diagnostic. 

In  the  last  few  years  there  has  been  extensive  Interest 
In  the  techniques  of  the  PTDS  and  PADS.  Since  the  Initial 
work  of  Davis  (18:515)  and  of  Boccara,  Fournier,  Amer  and 
co-workers  (19:130:  20:377;  21:519),  extensive  new 
applications  have  been  developed  by  Tam  and  co-workers 
(22:457:  23:436;  24:280;  25)  and  by  others  (26:381).  Qupta 
and  co-workers  have  demonstrated  the  usefulness  of  these 


techniques  for  combustion  diagnostics  (27:2663;  26:3673: 
29:303;  30:672;  31:4626).  The  two  techniques  are  closely 
related  to  each  other  and  satisfy  most  of  the  criteria 
enumerated  above  for  an  “ideal"  technique. 

Besides  the  widely  used  laser-Doppler  anemometry  (LDA) 
(32;  33:473)  and  laser-enhanced  ionization  (LEI)  (34:166) 
techniques  applied  to  the  flow  velocity  measurements  in  a 
flame,  other  techniques  such  as  PTDS  (35:532)  which 
overcomes  the  difficulty  associated  with  LDA  for  low 
velocities,  and  image  processing  (IP)  (36:1420)  have  been 
demonstrated . 

As  a  matter  of  fact,  no  ideal  technique  seems  to  exist 
and  all  of  the  techniques  have  relative  advantages  and 
disadvantages  over  each  other.  Table  I  shows  some  of  these 
characteristics  for  comparison  in  temperature  measurements 
of  flame  environment. 

The  extension  of  the  non-linear  techniques  to  high 
frequencies  (essential  to  an  understanding  of  turbulence 
phenomena)  has  suffered  from  the  lack  of  high  repetition- 
rate,  high  peak-power  laser  sources.  High-frequency 
thermometry  has.  to  date,  been  carried  out  using  fine-wire 
thermocouples,  Rayleigh  scattering  (although  both  methods 
are  extremely  difficult  to  apply  to  a  practical  flame 
system),  and  the  PTDS  technique  (37;  38). 

A  novel  technique  capable  of  high-frequency  thermometry 
along  a  line  is  the  thin  filament  pyrometry  (TFP)  technique 


TABLE  1 


Comparison  Between  Various  Optical 
Techniques  for  Thermometr/ 


Tech 

nlque 

Advantages 

Disadvantages 

OAKS 

High  spatial  resolution. 

High  sensitivity. 

High  peak-power  pulsed 
laser  needed. 

Expensive  Instrumenta¬ 
tion. 

Susceptibility  to  envi¬ 
ronmental  disturbances. 
No  applicable  to  mino¬ 
rity  species. 

LIB 

High  spatial  resolution. 
High  sensitivity. 

Simple  Instrumentation . 

Need  to  know  quenching 
rates 

BLIB' 

High  spatial  resolution. 

High  sensitivity. 

Lower  sensitivity  with 
molecules  with  long- 
lived  excited  states. 
High  laser  power  needed. 
Not  applicable  to  mole¬ 
cular  species  with  tran¬ 
sitions  In  the  UV  or 
violet  region. 

FAS 

2-D  resolution. 

Resolution  limited  by 
mlcrophon  bandwidth. 

Foor  sensitivity  using 
FZT  transducer. 
Interference  by  the 
acoustic  wave  passing 
through  the  flame  sheet. 

FADS 

High  2-D  spatial 
resolution. 

High  temporal 
resolution. 

Applicable  In 
sooting  tlames. 

Lower  sensitivity  than 
FTDS.  but  Increasing 
as  T  decrease. 

TABLE  1  (Cont'd) 

Comparison  Between  Various  Optical 
Techniques  lor  Thermometry 


Tech 

nlque 

Advantages 

Disadvantages 

PTDS 

High  2-0  spatial 
resolution. 

High  sensitivity. 
Applicable  In 
sooting  t lames. 

High  temporal 
resolution. 

that  has  been  recently  applied  successfully  by  Vllllmpoc  et 
al.  to  the  study  of  diffusion  and  premixed  flames  to  obtain 
time  resolved  temperature  profiles  (39;  40).  This  technique 
relies  upon  the  blackbody  emission  of  a  thin  ceramic 
filament  (SIC)  placed  In  a  combusting  flow  field.  The  major 
attributes  of  this  method  are  the  high-frequency  response, 
good  spatial  resolution,  and  the  ability  to  make  point 
temperature  measurements  along  the  line  of  the  cer6unlc 
filament  with  flames  In  the  range  1000‘K-2600‘K.  The 
technique  Is  limited  to  non-sootlng  flames  since  the 
emission  from  soot  can  Interfere  with  the  filament  emission. 

Finally,  the  flow  visualization  Is  an  essential  part  in 
the  combustion  study.  It  Is  used  (1)  to  Identify  the  scales 
of  processes  taking  place  In  different  regions  of  combusting 
flows.  (11)  as  an  aid  In  Interpreting  measurements  from  any 
diagnostic  technique,  and  (111)  to  provide  Insights  Into  the 
dynamics  of  the  mixing  and  transport  processes  Inside  of 
flames.  A  new  laser-sheet  lighting  technique  has  recently 
been  demonstrated  (41:159-174;  42;  43);  It  Involves  seeding 
the  central  fuel  jet  with  T1C14  vapor,  as  the  fuel  is  mixed 
with  air  the  TICI4  reacts  spontaneously  and  nearly 
Instantaneously  with  the  water  vapor  In  the  flame  to  form 
micro  size  TIO2  particles  and  HCl .  The  Mle  scattering  from 
the  T102  particles,  observed  at  right  angles  to  the  plane  of 
the  laser-sheet,  provides  a  planar  view  of  the  flow  field 
Inside  the  flame.  The  combination  of  two  perpendicular 


laser-sheetB.  aa  la  ahown  In  Fig.  1,  haa  given  two  bl- 
dlmenslonal  views  of  the  flame  structure  . 

Objective  and  Approach 

The  main  purpose  of  this  study  was  to  obtain  an 
experimental  characterization  on  temperature  at  high 
repetition  rates  for  a  CO-Hs  diffusion  flame  using  the  thin 
filament  pyrometry  technique.  The  present  measurement  has 
been  limited  to  the  near  field,  e.g.  the  region  within  20 
diameters  downstream  from  the  burner  exit. 

There  are  two  main  reasons  for  selecting  a  CO-H2 
diffusion  flame  In  this  study:  (1)  the  pyrolysis  of  high 
hydrocarbons  constituents  of  common  fuels  produces  CO  and 
Hz  which  react  with  the  ambient  air  In  a  fast  oxidation  In 
the  flame,  and  Cll)  using  CO  and  Hz  one  can  easily  control 
the  chemistry  of  the  flame  and  obtain  accurate  predictions 
for  other  fuels. 

Since  the  TFP  technique  has  been  used  In  a  Nz  diluted 
Hz  diffusion  flame  (39;  40)  this  study  was  undertaken  by 
using  a  Nz -co-Hz  mixing  as  a  previous  step  to  the 
thermometry  In  CO-Hz  mixing.  On  the  other  hand,  the  fact 
that  both  CO  and  Nz  have  the  same  molecular  weight  and 
practically  the  same  thermal  conductivity  and  viscosity 
coefficients  leads  one  to  expect  a  similar  flame  structure 
for  these  flames  Judging  from  the  fluid  mechanics  point  of 
view,  although  the  chemistry  Is  quite  different. 

Measurements  were  taken  at  various  axial  locations  of 
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B1 -dimensional  visualization  of  a  flame  In 
two  perpendicular  planes. 


the  flame  with  a  step  resolution  of  10  mm  using  different 
gas  mixtures  as  the  way  to  get  the  temporal -radial 
temperature  profiles  at  several  fluid  dynamics  regimes. 

Furthermore,  the  two  perpendicular  laser-sheets 
lighting  technique  was  used  to  visualize  the  flow  as  an  aid 
to  the  understanding  of  the  dynamic  flame  structure. 


II.  Ibeoo 


Thin  Filament  Pyrometry  Technique 

This  technique  relies  upon  the  blackbody  emission 
(radiation  resulting  from  the  thermally  excited  molecular 
oscillations  of  the  source)  of  a  ultrafine  0-SiC  (with 
excess  carbon)  filament  available  commercially  (44).  The 
filament  diameter  ranges  from  10  to  20  urn  (nominally  15  urn), 
and  the  thermal  conductivity  is  10  Kcal/m  hr“C  along  the 
filament  axis  at  room  temperature.  Because  of  its  small 
size,  the  convective  eat  transfer  coeficient  of  the 
filament  is  quite  large  and  it  can  respond  almost 
instantaneously  to  temperature  changes  in  its  surroundings; 
estimates  of  fiber  response  to  a  transient  heat  pulse  are 
between  1  ms  and  2  ms .  i.e.  •»>700Hz. 

The  thermal  conductivity  of  the  filament  is 
approximately  forty  times  lower  than  that  of  a  Pt-PtRh  wire 
which  drastically  reduces  the  spread  in  the  temperature 
profile  along  the  filament,  resulting  in  high  spatial 
resolution.  Because  the  emissivity  (ratio  of  the  excitance 
of  the  actual  source  at  each  temperature  to  the  ideal 
blackbody  excitance  at  the  same  temperature)  of  the  filament 
is  large  (0.88)  and  constant  with  respect  to  wavelength,  the 
fiber  can  be  considered  as  a  graybody  in  radiometric  terms, 
and  its  spectral  excitance.  Me,  (radiation  flux  emitted  per 
unit  area,  per  unit  wave  length  in  all  directions)  can  be 
directly  related  to  the  flame  temperature. 


Quantitative  measurements  with  the  filament  require 
that  its  temperature  response  be  calibrated  in  order  to 
establish  a  relationship  between  filament  emission  intensity 
and  adiabatic  temperature.  The  equations  describing  the 
blackbody  emission  from  such  a  filament  are  well 
established.  The  Planck  equation  which  describes  the 
spectral  distribution  of  excitance  for  an  ideal  blackbody 
at  each  wavelength,  X, 

2xhc*  1 

Me(X)  - - (1) 

X5  exp(hc/kXT)  -  1 

is  used  in  conjunction  with  the  Gouff6  equation  to  determine 
the  totai  spectral  flux  which  is  emitted  from  the  surface  of 
the  filament 

a  €o (Xo ,T)  Cl 

Me(Xo.T)  = -  (Watt/m)  (2) 

Eexp(C2 /XoT)  -  1] 

where  a  is  the  filament  surface  area,  €o  the  apparent 
emissivity  of  the  filament.  Cl  the  first  radiation  constant, 
C2  the  second  radiation  constant,  Xo  the  wavelength  in 
vacuum,  and  T  the  filament  temperature.  If  Eq.(l)  is  to  be 
used  for  quantitative  measurement  of  the  temperature  of  the 
filament,  the  emissivity  of  the  filament  must  be  known  as  a 
function  of  temperature.  The  conventional  approach  to 


calibrating  the  emissivity  of  a  material  is  to  heat  the 
given  material  to  a  known  temperature,  measure  its 
excitance,  and  compare  it  to  that  of  a  blackbody  at  the  same 
temperature,  but  due  to  the  small  size  filament  employed, 
this  approach  is  impractical.  Therefore,  as  established  by 
Goss  et  al.  (45)  an  alternative  approach  of  measuring  the 
reflectance  (ratio  of  the  power  reflected  to  the  power 
incident)  of  the  filament  as  a  function  of  temperature  can 
be  used. 

The  reflectance  can  be  related  to  the  directional 
emissivity,  G(0,S),  by  Kirkoff’s  relationship 

6(0,*)  =  1  -  p(0,4)  (3) 

where  P(0,*)  is  the  directional  hemispherical  reflectance 
for  the  spherical  coordinates  0  and  4.  In  general,  the 
required  accuracy  in  reflectance  for  a  desired  accuracy  in 
emissivity  accuracy  is  given  by  the  following  relation 


AP  e  4€ 


where  Ap  and  d€  represent  the  corresponding  errors.  For  a 
highly  emissive  sample  such  as  the  SiC  filament,  with 
emissivity  of  0.88,  the  error  in  the  emissivity  measurement 


based  on  a  1%  error  in  the  reflectance  is  0.13%.  Thus, 


reflectance  measurements  (if  carried  out  correctly,  i.e. 
measuring  the  bi-directional  reflectance  distribution 
function  at  all  collection  angles)  can  be  very  accurate 
measure  of  emissivity.  The  concern  with  the  SiC  fiber  is 
not  to  measure  the  absolute  value  of  the  emissivity  at  all 
temperatures  but  rather  to  establish  whether  the  emissivity 
remains  constant  in  the  range  of  interest.  This  requires 
that  the  change  in  the  emissivity  (not  the  absolute  value) 
be  monitored  as  a  function  of  temperature.  This  can  be 
accomplished  with  the  experimental  setup  shown  in  Fig.  2. 

The  chopped  output  from  a  HeNe  laser  is  monitored  by  a 
photomultiplier  tube,  PMT,  and  lock-in  amplifier.  The 
bundle  of  filaments  are  heated  by  a  flame  to  a  maximum 
temperature  of  ~2200‘’K.  Since  the  filament  is  a  graybody 
and  since  no  measurable  change  is  observed  in  the  emissivity 
at  the  HeNe  wavelengths,  it  can  be  concluded  that  the 
emissivity  of  the  filament  indeed  remains  constant  over  the 
temperature  range  of  interest.  This  implies  that  the 
excitance  of  the  filament  can  be  predicted  by  the  blackbody 
equations.  In  practice,  however,  a  detector  of  limited 
spectral  bandwidth  is  used  to  detect  the  excitance  and, 
thus,  the  blackbody  results  must  be  convolved  over  the 
detector  response  curves  to  predict  the  measured  response  of 
the  filament. 

The  equation  (2)  must  be  multiplied  by  the  detector 
response  and  integrated  over  this  region.  Thus,  the 


detected  signal  is  given  by 


lex  p 


kex  p 


'A2 

Me(A.T)  R(X)  Tr(X)  dA 
Ai 


(5) 


where  Me(A,T)  is  the  graybody  function,  R(A)  the  response 
function  of  of  the  InGaAs  photodiode,  Tr(A)  the  transmission 
characteristics  of  the  optics  employed  in  the  experiment, 
and  kexp  an  experimental  constant  which  takes  into  account 
the  collection  efficiency  of  the  optical  system. 

To  simplify  the  calibration  procedure,  the  emission 
from  the  filament  at  a  known  temperature  is  measured,  and 
all  measurements  are  referenced  to  this  value.  This  allows 
experimental  factors  (alignment,  collection  efficiency, 
fiber  diameter,  etc. )  to  be  normalized.  Moreover,  the 
normalization  between  the  different  filaments  used  through 
the  same  experiment  is  provided  by  ratioing  the  measurements 
at  an  axial  location  in  the  flame  where  the  emission  from 
the  filament  does  not  change  in  time. 

Finally,  the  intensity  data  are  converted  to 
temperature  using  a  calibration  curve.  This  curve,  shown  in 
Fig.  3  (45).  is  calculated  by  convolving  the  detector 
response  over  the  blackbody  excitance;  the  curve  is  fit  by  a 
mathematical  power  series  expansion  of  intensity  to 
determine  the  temperature  and  normalized  to  the  emission  at 


a  known  temperature  and  axial  location. 


Because  the  temperature  of  the  filament  may  be 
different  from  that  of  the  surrounding  gas  due  to  radiation 
losses,  a  detailed  heat-balance  analysis  must  be  conducted 
for  accurate  prediction  of  the  gas  temperature.  The  results 
of  such  an  analysis  for  the  flame  conditions  reported  by 
Goss  et  al.  (45)  indicated  that  the  radiation  correction 
varies  from  20“K  at  1000“K  to  SOO^K  at  2600‘’K. 

As  can  be  seen  from  Fig.  3,  the  detected  emission  from 
the  filament  varies  greatly  with  temperature.  Because  of 
this  variance  along  with  the  limited  dynamic  range  of  the 
digitization  system  (255  to  1),  the  temperature  uncertainty 
can  be  large  (especially  at  low  temperature).  Assuming  a  1% 
uncertainty  in  the  intensity  measurements,  the  temperature 
uncertainty  at  1100  “K  is  »80  *K  and  at  2300  "K  is  -8  "K. 
LaserrSheet  Lighting  Technique 

The  laser-sheet  lighting  technique  involves  seeding  the 
central  fuel  jet  with  TiCl4  vapor.  As  the  fuel  is  mixed 
with  the  air,  the  TiCl4  reacts  with  the  water  vapor  formed 
in  the  combustion  process  to  form  Ti02  particles  and  HCl 
according  with  the  stoichiometry 

TiCl4  +  2H20  --->Ti02  +  4HC1  (6) 

The  Mie  scattering  from  the  Ti02  particles  is  observed 
as  a  green  color  at  right  angles  to  the  planes  of  the  green 
Nd:YAG  laser-sheets  which  are  passed  vertically  through  the 


center  of  the  flame  and  horizontally  at  the  desired  axial 
location.  The  technique  Is  simple  and  easy  to  Implement 
requiring  only  a  laser  and  camera. 

The  Introduction  of  TICI4  to  the  flow  does  not  disturb 
the  combustion  process  since  the  reaction  of  TICI4  vapor  and 
H2O  vapor  occurs  nearly  Isothermal ly  and  the  amount  of  TICI4 
added  to  both  the  dry  air  and  the  dry  central  fuel  Is 
negligible. 


III.  Experimental  Apparatus 


Experimeritai  Setup  for  the  TFP  Technique 

The  optical  arrangement  employed  in  this  technique  is 
shown  in  Fig.  5.  An  axisymmetric  fuel  nozzle  of  10  mm 
diameter  was  mounted  at  the  center  of  a  co-axial  air  jet  of 
2b4  mm  diameter;  the  low  velocity  (16  SCFM.  cm/s)  air  is 
introduced  through  the  outer  tube  about  the  fuel  tube  to 
reduce  the  effect  of  room  air  disturbances.  The  burner  was 
mcjunted  on  a  motor-driven  translation  stage  to  permit 
vertical  movements.  A  brief  description  of  the  main 
instruments  employed  in  this  technique  and  their  principal 
characteristics  follows. 

Rotating  mirrop 

In  the  thin  filament  pyrometry  technique  the  light  from 
the  filament  is  imaged  by  a  44  cm  focal  length  lens  onto  a 
rotating  mirror,  which  sweeps  the  image  of  the  filament 
across  a  single  detector  through  a  50  cm  focal  length  lens. 
This  mirror-detector  combination  converts  a  spatial  scan 
into  a  temporal  one. 

The  Lincoln  Laser  Corp. .  model  S-225-030-XLIM. 
rotating  mirror  scanner  (RMS)  is  actually  a  wheel  consisting 
of  ten  different  faces  which  rotates  at  speeds  up  to  bOO  Hz. 
The  wheel  in  these  experiments  rotated  at  a  100  Hz  rate 
driven  by  a  variable  frequency  motor  control  Lincoln  Laser 
Corp.  model  VFC-2;  with  the  ten  faces,  the  image  was  swept 
by  the  detector  at  an  effective  rate  of  1000  Hz. 


The  image  from  the  rotating  mirror  is  swept  across  a 
In-Ga-As  photodiode.  RCA  model  C30980E,  whose  spectral 
response  is  shown  in  pg.  19,  and  whose  principal  technical 
specifications  are: 

-quantum  efficiency  at  1300  nm  and  1550  nm,  70%  typ. 

-total  dark  current  (Id),  150  nA  typically. 

-noise  current  (In)  at  f  =  10  kHz,  0.3  pA/Hz%  typ. 

Due  to  the  spectral  response  which  is  limited  to  the 
near  IR,  only  temperatures  in  excess  of  1100*K  could  be 
measured.  However,  with  a  suitable  choice  of  IR  detector 
and  optics,  this  lower  temperature  limit  could  be  extended. 
Digitizer 

A  Lecroy  TR8818  Transient  Recorder  is  used  to  digitize 
the  analog  signal  from  the  detector  and  store  it  in  MM8103 
ECL  memory  units  (32  K  each  one)  for  subsequent  transmission 
through  a  custom  built  interface  to  a  ModComp  7870  Classic 
computer,  for  storage  and  data  analysis  (in  which  the  8  bit 
data  from  the  LeCroy  are  unpacked  into  16  bit  data  for  the 
ModComp).  The  synchronization  pulse  from  the  PIN-diode  is 
delayed  and  used  to  open  a  200  us  window  centered  on  the 
sweep  of  the  filament  image;  this  delay  allows  the  image 
signal  to  be  located  at  the  optimum  point  inside  the  200  us 
window.  The  image  of  a  6  cm  length  along  the  filament  is 
recorded  in  512  discrete  samples  during  the  200  us  window  at 
the  rate  of  2.56  MHz.  giving  an  effective  spatial  resolution 


of  120  urn  per  sample  point.  A  pulse  generator  (General 
Radio  model  1340),  was  used  to  provide  the  512  pulses  which 
carry  at  the  leading  edge  the  signal  from  the  In-Ga-As 
detector . 

The  model  TR8818  digitizer  is  optimized  for  digitizing 
and  recording  high  frequency  signal,  evidenced  by  its  wide 
bandwidth  (DC  to  >100  MHz),  fast  sampling  rate  and  its  high 
dynamic  accuracy  even  for  signal  frequencies  approaching  the 
Nyquist  limit.  The  signal  to  be  digitized  is  connected  via 
an  amplifier  ( Lecroy  6102)  for  gain  control.  Once  armed  by 
programmed  command  through  computer  keyboard,  the  unit 
digitizes  the  input  signal,  with  8-bit  resolution  and  up  to 
200  MS/s  digitization  rate,  and  loads  its  memory. 

Triggering  Lasej; 

A  HeNe  laser  (Spectra  Physics,  model  145-01)  and  a  PIN- 
photodiode  providing  the  necessary  synchronization  signals 
were  used  for  proper  scanning  and  digitization  of  the  image. 
The  laser  was  used  with  the  power  supply  Spectra  Physics 
model  248. 

Experimental  Setup  for  the  LrS  Lighting  Technique 

The  experimental  arrangement  for  this  technique  is 
shown  in  Fig.  6.  A  pulsed  frequency  doubled  Nd:YAG  laser 
(Quanta-Ray  DCR)  was  used  as  the  laser  source  to  produce  the 
horizontal  and  vertical  laser  beam  sheets.  The  principal 
characteristics  of  this  laser  are: 


-pulse  width.  10  ns. 


-pulse  energy,  70  mJ/pulse  with  amplifier  oscillator, 
-optimum  pulse  repetition  rate,  10  pulses/s. 

Although  the  power  of  the  Nd:YAG  is  adjustable,  the 
higher  range  of  operation  was  used  to  provide  better 
illumination  to  the  internal  structure  of  the  flame  and 
obtain  the  best  definition  on  the  photo. 

A  Hasselblad  500  camera  was  used  whose  shutter  was 
triggered  by  the  laser  firing  using  an  exposure  time  of 
1/500  s  and  a  diaphragm  aperture,  fit,  of  4  and  5.6. 

In  this  experiment  both  the  laser’s  oscillator  and 
amplifier  were  used  to  generate  the  required  output  energy. 
A  burner  with  similar  geometry  to  the  one  used  in  the  TFP 
technique  was  employed  in  this  experiment. 

Since  HCl  is  formed  in  the  TiCl4  reaction  with  the 
water  vapor,  some  precautions  must  be  taken  when  applying 
the  technique.  Since  an  extremely  small  amount  of  TiCl4  is 
used  to  seed  the  flow,  the  concentration  of  the  HCl  formed 
is  very  small.  However,  care  must  be  taken  to  vent  the 
exhaust  products  from  the  laboratory. 


IV.  Experimental  Procedure 


Experimental  Procedure  for  the  TFP  Technique 

In  the  thin  filament  pyrometry  technique,  using  the 
optical  arrangement  shown  in  pg.  23,  radial  temperature 
profiles  of  the  flame  as  a  function  of  time  were  made  at 
different  axial  locations,  from  the  nozzle  exit  to  200  mm 
above  the  nozzle  in  10  mm  increments.  The  test  conditions 
are  summarized  in  TABLE  II. 

The  filament  was  suspended  on  the  centerline  of  the 
diffusion  jet  by  a  holder  which  supplies  just  enough  tension 
to  keep  the  filament  from  mechanical  oscillations  (probable 
at  relatively  high  jet  velocities)  but  not  too  much  as  to 
decrease  drastically  its  survivability  at  high  temperatures 
where  its  tensile  strength  is  reduced.  The  emission  from 
the  filament  was  swept  across  the  photodetector  by  the 
spinning  mirror  and  sent  to  the  computer  after  the  signal  is 
digitized.  The  digitiser  used  in  the  study  stores  384 
images  digitizing  each  image  in  512  records.  Depending  on 
the  intensity  of  the  emission  three  different  electronic 
ranges  (50  mV,  100  mV  and  250  mV  full  scale)  were  employed 
for  digitization;  in  all  ranges  the  full  scale  corresponded 
to  255  counts. 

Cal i bra t ion  Procedure 

The  environme'it  in  which  the  filament  is  placed  is 
harsh.  Each  filament  could  last  only  a  few  minutes 
(depending  on  axial  location.  Re  number  and  gas  mixture 


TABLE  M 


Test  Conditions  for  Temperature  Measurementj 
of  H2 -CO-N2  and  H2 -CO  Fiames 


Test  Conditions 

Experiment 

B 

C 

D 

E 

F 

Fuel 

H2  - 
-CO-N2 

H2- 

-CO-N2 

H2- 

-CO-N2 

H2- 

-CO-N2 

H2  - 

-CO-N2 

H2  -CO 

Re* 

1717 

2099 

854 

857 

2264 

857 

Exit  velocity 

V  (m/s) 

2.69 

3.92 

1.60 

2.87 

4.90 

2.87 

H2 /(CO+N2 ) 


H2  /CO 


0.029  0.497  0.504  3.391  0.998 


3.391 


(  •♦  )  Re  -  40/ndl/  ;  0  =  volumetric  flow  rate  in  cm3/s. 
d  -  burner  diameter  in  cm.  v  -  kinematic  viscosity  of  fuel 
gas  at  298  “K  in  cm2 /s . 

ft*)  R.i.p.m.  is  based  on  the  rotameter  reading  in  liters 
per  minute  calibrated  for  C2H2 . 


employed),  so  several  fibers  were  used  in  each  experiment. 
The  correlation  between  measurements  made  with  different 
fibers  in  each  experiment  was  accomplished  by  measuring  at 
10  mm  above  the  nozzle  with  all  filaments.  Calibration 
measurements  at  40  mm  downstream  were  made  for  each  matrix 
using  pure  H2  flame.  The  flame  background  excitance 
was  measured  at  20  mm  steps  for  correcting  the  filament 
measurements.  The  goal  of  this  process  is  to  obtain  the 
excitance  at  the  adiabatic  temperature  of  the  flame. 

Data  Analysis 

Data  analysis  of  the  experimental  data  was  carried  out 
on  the  Modcomp  Classic  computer.  This  analysis  required 
several  steps  which  included: 

-noise  reduction  by  averaging  six  points  when  the 
difference  between  two  consecutive  records  in  each  image  is 
bigger  than  5  counts. 

-correlation  of  all  measurements  made  with  the 
different  fibers,  by  ratioing  the  maximum  values  obtained  at 
10  mm  above  nozzle  with  the  first  to  the  subsequent 
filaments . 

-conversion  between  the  different  electronic  scales 
used  in  the  digitizer  by  converting  the  number  of  counts  of 
the  digital  signal  to  actual  milivolts. 

-background  compensation  of  the  signal  by  subtraction 
of  averaged  background  from  each  record  through  all  images. 

-calibration  with  pure  H2  flame  by  ratioing  each  record 


in  each  image  to  the  maximum  value  obtained  at  40  mm  above 


nozzle,  and  converting  the  ratioed  intensity  to  temperature 
using  the  fitted  calibration  curve  relative  to  the 
stoichiometric  temperature  of  a  Hz  flame  convolved  with  the 
detector  response  (45). 

This  process  required  a  modif icationt  of  available 
computer  programs  to  include  the  different  factors  involved. 
The  actual  temperature,  TEMP,  was  obtained  from  the 
intensity,  INT,  by  using  the  calibration  curve  and  assigning 
the  limits 

TEMP  =  2600  “K  when  INT  =  1 

TEMP  =  1100  “K  when  INT  ^  0  (7) 

where 

Countsi  Fwi  (mV/counts  )i 

INT  = - 

Max  Countscal  Fwcai  (mV/counts )c a i 


Bkg  Countsi  (mV /counts )Bkgi 


Max  Countscai  (mV/counts)cai 

Countsi /(counts/mV )i  (max  countswo/max  countswi  ) 


Max  Countscai /(counts/mV)cai  (max  countswo /max  countswcai  ) 


Bkg  Countsi /( counts/mV )Bkgi 


Max  Countscai /(counts/mV)cai 


[(Countsi  /  (counts/mV )i )/max  countswi  ] 


[(Max  Countscai /(counts/mV )c a  1 /max  countswcai  j 


Bkg  Countsi /( counts/mV )Bkgi 

-  (8) 

Max  Countscai /(counts/mV )cai 


where 

Countsi 

Max  Countscai 
( mV/counts  )i 

(mV/counts )cai 

(mV/counts  )Bkgi 

Fwi 


Fwc  a  1 


counts  at  location  i  in  each  record  of 
every  image . 

maximum  counts  at  calibration  measurement, 
milivolts  per  count  of  digitizer  scale  at 
location  i. 

milivolts  per  count  of  digitizer  scale  in 
calibration  measurement. 

milivolts  per  count  of  digitizer  scale  at 
location  i  in  background  measurement, 
conversion  factor  of  fiber  at  location  i 
to  location  at  0  mm  (ratio  of  maximum 
counts  at  10  mm  above  nozzle  with  the 
fibers  used  at  locations  0  mm  to  i  mm), 
conversion  factor  of  fiber  in  calibration 
measurement  to  location  at  0  mm  ( ratio  of 
maximum  counts  at  10  mm  above  nozzle  with 
the  fiber  used  at  location  0  mm  to 
calibration  measurement  fiber) . 


Bkg  Countsi 


counts  at  location  i  in  each  record  of 


every  Image  in  background  meaaurement. 

As  has  been  pointed  out  by  Goss  (45)  the  measured 
filament  temperature  has  to  be  corrected  by  radiation  losses 
to  obtain  the  surrounding  gas  temperature.  As  the 
temperature  Increases  the  blackbody  emission  Increases  as 
T4  (according  with  the  Stefan-Bolztmann  formula)  while  the 
convective  heat  transfer  Increases  only  as  T.  thus  the 
correction  term  becomes  large  as  the  temperature  Increases . 
and  Its  magnitude  depends  on  the  flow  conditions  and 
constituents  as  well  as  the  temperature  Itself  (through  the 
influence  of  temperature  on  the  heat  transfer  coefficient). 
So.  the  correction  term  has  to  be  evaluated  each  time  the 
flow  conditions  or  fuel  types  are  changed.  This  correction 
factor  was  not  introduced  in  the  three-dimensional  plots. 

A  graphic  printer.  Tektronix  model  4012.  was  used  to 
plot  the  temporal  variation  of  the  temperature  across  the 
radial  direction  of  the  flame,  and  the  Intensity  contours  In 
a  tlme-versus-radlal  distance  plot. 

Because  of  the  nature  of  the  plots  obtained,  the 
excltance  data  were  spectrally  analyzed  by  using  a  Fourier 
Transform  analysis  computer  program,  to  show  the 
distribution  of  the  temperature  variation  In  the  frequency 
domain  along  the  radial  direction  In  the  flame. 


V.  Results  and  Discussion 


f'iow  Visualization  Studies 

Yule  et  al.  have  made  detailed  studies  of  the  structure 
of  a  transitional/turbulent  propane  Jet  flame  with  a  burner 
setup  very  similar  to  the  present  studies  (46).  They 
identified  two  regions  containing  organized  structures.  One 
region  located  outside  of  the  flame  surface,  consisted  of 
large  slow  moving  vortices  that  occurred  at  a  frequency  of 
about  30  Hz.  They  referred  to  this  zone  as  the  outer 
preheat  zone.  Another  zone,  referred  to  as  the  inner 
preheat  zone,  was  located  inside  the  flame  at  the  potential 
core  interface.  These  types  of  structures  shown  in  Fig.  7a 
&  7b.  are  clearly  identified  in  the  visualizations  of  the 
H2 -CO  flame  investigated  in  this  study.  The  structures 
inside  the  visible  flame  surface  mark  the  fuel  and  HzO 
combustion  product  interface,  and  the  large  structures 
outside  the  flame  surface  mark  the  air  and  H2O  product 
interface . 

Studies  by  Roquemore  et  al .  (47)  have  shown  that  at  low 
and  intermediate  velocities,  the  toroidal  vortices  observed 
outside  the  visible  flame  surface  play  a  dominating  role  in 
determining  the  flame  structure.  The  rotation  of  the 
toroidal  vortices  can  pull  the  flame  radially  outward 
forming  a  flame  bulge.  The  frequency  of  the  outer  vortices 
corresponds  to  the  flame  flicker  frequency  of  •»12  Hz.  and  is 
relatively  independent  of  the  fuel  type  and  central  nozzle 
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design  for  a  wide  range  of  fuel  velocities.  The  outer 
structures  are  believed  to  be  established  by  a  buoyancy 
driven  shear  layer. 

The  turbulent  characteristics  of  the  jet  diffusion 
flame  are  most  clearly  demonstrated  by  the  fuel/water 
product  interface  occurring  inside  the  flame  surface.  An 
initial  instability  wave  forms  in  the  potential  core  region 
of  the  jet  at  very  low  fuel  flow  rates.  This  inner 
instability  can  develop  into  coherent  vortices  as  the  fuel 
velocity  increases.  The  flame  spreading  angle  increases  at 
the  axial  location  where  the  vortices  start  to  coalesce  much 
like  a  non-reacting  shear  layer  grows  when  vortices 
coalesce.  At  the  low  flow  rates  of  the  flame  investigated 
in  this  study,  the  inner  vortices  are  not  well  form  as 
illustrated  in  Fig.  7a  &  7b. 

The  sheet  lighting  technique  illustrates  the  internal 
and  external  flame  structures  but  does  not  show  the 
temperature  field.  The  TFP  technique  is  well  suited  for 
this  purpose  as  discussed  below. 

TFP  Technique 

The  development  of  the  present  experiments  has 
indicated  the  necessity  for  enhancing  the  performance  of  the 
filament.  The  low  survivavility  of  the  filament  under  the 
conditions  used  here  -in  contrast  with  the  reported  results 
using  H2 -N2  diffusion  flame  (38;  39)-  could  be  due  to  some 
combined  chemical  effects  of  the  CO  and  O2  over  the  SiC 


fiber.  This  situation  leads  to  a  lower  reliability  of  the 
results  and  make  its  applicability  to  in  situ  measurements 
uncertain.  A  recently  developed  ceramic  fiber  (boron 
nitride  coated  by  Si02  and  Si3N4)  seems  to  be  more  resistant 
to  the  oxygen  action  increasing  its  lifetime  in  flame 
environment . 

The  addition  of  extended  memory  to  the  digitizer  and 
computer  will  allow  an  expanded  frequency  analysis  to  high 
frequencies,  and  make  enough  point  and  line  measurements 
along  the  time  to  get  averaged  and  probability  distribution 
function  (PDF)  data  as  a  spatial  and  temporal  function. 
Temperature  Profiles 

The  3-D  temporal  distribution  of  temperature  profiles 
along  the  radial  direction  taken  at  various  axial  locations 
are  shown  from  Fig.  8  to  Fig.  18  corresponding  to  the 
experiment  F  (pure  CO-H2  flame).  Included  in  these  figures 
are  the  respective  iso-intensity  contour  lines  in  a  time 
versus  radial  distance  graph.  The  reader  should  keep  in 
mind  that  the  displayed  temperature  fields  represent  the 
passage  of  the  flame  structure  through  a  single  axial 
location  as  a  function  of  time,  so  the  y-axis  in  these  plots 
is  time  and  not  space  as  in  the  photograph  of  the  flame. 

The  evolution  of  a  large  scale  buoyancy  driven 
structures  (flame  bulges)  can  be  traced  from  these  plots. 
Just  above  the  nozzle  exit  the  cold  core  of  unburnt  fuel  is 
surrounded  by  a  hot  reacting  flame  sheet  (surface)  resulting 
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Fig.  8  Temperature  Profiles  of  CO-Hz  Diffusion 
Flame  and  Intensity  Contour  Lines  at 
Axial  Location. 


from  the  combustion,  at  the  stoichiometry  rate,  of  fuei  and 
the  co-flowing  air.  The  reaction  zone  is  well  defined  and 
steady  due  to  the  molecular-transport  nature  of  the  mixing 
processes  in  this  region.  The  flame  structure  display  no 
indication  of  any  low  frequency  oscillation  but  a 
comparatively  high-frequency  oscillating  structure  starts  to 
be  developed  in  the  internal  core  affecting  the  external 
flame  layer. 

At  20  mm  axial  location,  the  flame  has  increased  in 
size  radially  as  well  as  the  peak  temperature.  A  low- 
frequency  sinusoidal  oscillation  of  the  outer  flame  boundary 
is  occurring,  and  the  internal  high  frequency  is  fairly  well 
developed;  this  high  frequency  seems  to  be  in  the  same  range 
of  that  reported  by  Chen  (48)  for  the  first  sub-harmonic  in 
the  near  field  of  cold  round  jets.  The  reaction  sheet  has 
doubled  its  width  and  acts  as  a  damping  wall  of  the 
frequency.  Notice  that  the  temperature  gradient  is  sharper 
at  the  external  side  of  the  reaction  sheet  than  at  the 
internal  one.  The  increase  in  the  hot  zone  size  is  believed 
to  be  due  to  the  high  diffusivity  of  the  H2  molecule.  The 
contour  of  the  fuel  jet  was  constructed  to  give  a  flat  top 
velocity  profile  at  the  nozzle  exit  producing  a  very  thin 
shear  layer  inside  the  flame  surface;  this  is  the  origin  of 
local  instabilities  that  travel,  becoming  small  vortices. 
These  vortices  grow  and  coalesce  downstream  producing  a 
mixing  process  between  the  hot  reacting  layer  and  the  cold 


fuel  core.  These  vortices  together  with  the  high  viscosity 
of  this  flame  reaction  zone  can  be  responsible  for  damping 
oscillating  structure  and  the  smoother  temperature  gradient 
outside  tVie  elastic  flame  surface. 

At  40  mm  axial  location  the  jet  has  radially  expanded 
to  more  than  twice  its  initial  value  and  the  external 
oscillations  have  increased  markedly  in  amplitude  changing 
to  a  bell-like  shape.  Because  of  the  difference  in  the 
velocities  of  the  fuel  and  air  jets,  the  outside  shear  layer 
tends  to  transport  the  air  into  the  fuel  jet  promoting 
mixing.  The  reaction  zone  enlarges  as  more  and  more  air  is 
mixed  with  the  fuel  and  undergoes  combustion,  resulting  in  a 
rise  in  the  temperature  with  increasing  the  height  from  the 
nozzle  exit.  This  expansion  process  is  continued  at  60  mm 
and  80  mm.  changing  to  a  butterfly- like  shape.  The  flame 
surface  wraps  out  over  itself  which  results  in  the  formation 
of  a  mushroom-like  structure  that  engulfs  the  external  air 
and  forms  a  buoyancy  driven  toroidal  vortex.  This  bulge 
will  be  reproduced  at  16-17  nozzle  diameters  (D)  downstream, 
as  can  be  seen  on  Fig.  7  (obtained  using  the  2  laser-sheets 
lighting  technique).  At  this  bulge  location  the  diameter 
reaches  its  maximum  width  of  ^4  D. 

At  the  100  mm  axial  location  the  flame  speeds  up  in  the 
axial  direction  becoming  thinner  where  the  cold  central  jet 
is  now  unobservable.  The  mixing  process  drives  an  almost 
uniform  increasing-decreasing  process  of  the  radial  size. 


with  the  highest  temperature  occurring  at  the  maximum  width 
point.  This  process  continues,  driving  another  vortex  which 
feeds  the  flame  surface  with  external  air,  becoming  a  new 
bulge  approximately  at  160  mm  axial  location.  This 
interplay  between  the  bulges  and  the  coexisting  vortex 
structures  is  responsible  for  the  low  frequency  oscillations 
of  the  flame. 

Above  180  mm  axial  location  the  room  disturbances  start 
to  affect  the  flame  structure  giving  it  a  three-dimensional 
character  (flickering)  with  the  peak  temperature  decreasing 
once  the  majority  of  fuel  has  been  consumed  and  mixed. 

Higher  axial  temperature  profiles  were  not  attempted. 

The  height  at  which  the  toroidal  vortices  are  formed  is 
a  strong  function  of  the  co-annular  air  velocity,  being 
almost  unaffected  by  the  fuel  composition  (49).  Rough 
measurements  have  indicated  25-30  Hz  for  this  low  frequency 
oscillation,  while  the  high  frequency  oscillations  seems  to 
be  in  the  310-325  Hz  range. 

The  graphs  shown  in  Appendix  A  (from  Fig.  A-1  to  Fig. 
A-11)  corresponding  to  the  experiment  D,  reveals  a  similar 
behavior  when  a  H2-CO:N2  diffusion  flame  is  used.  All  of 
this  means  that  these  phenomena  are  predominantly  driven  by 
the  fluid  dynamic,  and  the  kinetic  reaction  is  initially 
confined  to  the  flame  sheet. 

The  Appendices  B,  C  and  D  include  the  graphs 
corresponding  to  the  experiments  E,  B  and  A  respectively. 


All  of  them  show  the  high-frequency  osci  i  1  at  iotis  in  the 
internal  fuel  jet. 

In  the  experiment  E  for  which  the  exit  velocity  is  4.9 
m/s  (notice  that  in  experiment  F  this  velocity  was  2.87 
m/s  )  .  the  external  side  of  the  shear  layer  starts  to  develop 
the  low-frequency  oscillations  at  40  mm  axial  location 
instead  of  20  mm  as  in  experiment  F.  The  cold  central  jet 
is  still  observable  at  200  mm  and  the  temperature  gradient 
is  roughly  equal  at  both  sides  of  the  shear  layer.  The 
general  flame  structure  shows  a  similar  behavior  to  the 
experiment  F. 

In  the  experiment  B  which  exit  velocity  is  3.92  m/s  the 
reaction  zone  is  well  defined  and  steady,  and  no  indication 
of  low  frequency  oscillations  is  displayed  up  to  60  mm  where 
the  radial  size  has  increased  twice.  At  70  mm  above  nozzle, 
the  butterfly-like  structure  indicates  that  the  bulging 
process  is  being  developed  in  the  flame. 

In  order  to  check  that  the  high-frequency  oscillations 
were  not  due  to  the  experimental  fixture,  a  chamber  (used  as 
pressure  wave  filter)  was  settled  to  the  fuel  feeding  lines 
using  the  gas  mixture  of  experiment  F.  Fig.  19  clearly 
shows  that  these  oscillations  are  an  internal  characteristic 
of  the  flame. 

As  a  way  to  investigate  if  higher  frequency 
oscillations  are  present  at  the  internal  structure  of  the 
flame  the  experiment  B  was  performed  using  higher  scanning 


Fig.  19 


Detail  of  Intensity  Contour  Lines  of 
CO-H2  Diffusion  Flame  at  40  nun 
Axial  Location. 


rates,  Same  results  were  obtained  at  different  effective 
swepts  up  to  0000  Hz. 

Analysis 

To  determine  the  frequencies  due  to  the  oscillating 
structure  obtained  from  the  intensity  contour  lines,  a 
Fourier  Transform  analysis  of  t,he  data  was  carried  out.  The 
3-D  plots  of  the  oscillation  spectrum  at  various  radial  and 
axial  locations  of  the  flame  are  shown  from  Fig.  20  to  Fig 
24  corresponding  to  the  experiment  F  (pure  CO-H2  flame).  As 
can  be  seen,  the  oscillations  when  developed  are  maintained 
along  the  axial  direction. 

Because  of  the  limitation  of  the  memory  load  cf 
digitizer  and  computer  used  in  the  analysis,  only 
frequencies  up  to  128  Hz  could  be  determined.  For  that 
reason,  it  could  not  be  established  the  exact  meaning  of  the 
plotted  frequency  set  around  84  Hz,  and  the  graphs  indicate 
only  tiiat  higher  frequencies  occur. 

At  30  mm  downstream  the  typical  low  frequency 
associated  with  the  buoyancy-driven  structures  starts  to 
develop;  at  90  mm  and  160  mm  axial  locations  the  low 
frequency  is  split  into  two  parts  radially  shifted.  This 
can  be  associated  to  the  bulges  developed  at  those  heights; 
the  surface  wraps  out  over  itself  forming  two  concentric 
flame  .sheets  during  a  limited  axial  range  until  tlie  bulge 
bec<>me.s  a  separate  portion  of  the  flame,  and  it  speeds  up 
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Frequency  Distribution  Profiles  of  CO-H2 
Diffusion  Flame  at  40  nun  and  60  mm  Axial 
Location. 


Fig.  25  shows  a  more  complicated  frequency  behavior  for 


VI.  CgnslysisQS  SeSfiOHDSDdSl'lfiDS 


The  thin  filament  pyrometry  technique  has  been 
demonstrated  to  be  capable  of  measuring  the  temporal-spatial 
temperature  distributions  in  a  non-sooting  CO-H2  diffusion 
flame  at  high  repetition  rates  up  to  5000  Hz.  Its  high 
spatial  and  temporal  resolution  and  its  capability  to  make 
point  temperature  measurements  along  a  line  allows  this 
technique  to  be  a  powerful  tool  for  a  quick  qualitative 
characterization  of  the  flame  identifying  its  critical  axial 
locations.  In  addition,  it  can  provide  a  quantitative 
approach  to  the  spatial  and  temporal  temperature 
distributions  as  well  as  to  the  frequency  spectrum  of  the 
fluid  dynamic  processes  related  with  the  flame. 

New  coated  boron  nitride  ceramic  fibers  will  increase 
its  survivability  under  the  severe  combustion  environment 
leading  to  a  higher  fiability  in  the  quantitative 
measurements.  A  suitable  choice  of  an  IR  detector  and 
optics  will  extend  the  lower  temperature  measurement  limit, 
set  at  1100  'K  at  present. 

Future  improvements  including  radiation  losses 
correction,  and  digitizer  and  computer  loading  memory 
increase  will  allow  a  higher  precision  in  the  temperature 
measurement,  a  higher  scanning  rate  with  the  consequential 
wider  frequency  analysis,  and  the  acquisition  of  sufficient 
amount  of  data  to  apply  a  statistical  analysis  and  determine 
average  values  with  its  PDF  along  the  time,  radial,  and 


axial  locations. 


An  intensive  theoretical  analysis  of  these  results  must 
be  carried  out.  Questions  like  why  the  flame  is  not 
stoichiometric  and  what  is  its  correlation  with  the  local 
diffusion  heat  losses;  where  is  the  largest  gradient  of 
temperature  and  why,  and  its  correlation  with  the 
temperature  peaks;  where  are  these  located  in  the  radial 
direction;  what  are  the  local  changes  of  the  thermal 
diffusion  coefficient;  why  is  the  rate  of  temperature  change 
different  at  the  internal  and  external  sides  of  the  reaction 
layer,  etc.,  need  to  be  answered.  The  TFP  technique  needs 
to  be  increased  in  its  quantitative  precision  through  a 
detailed  theoretical  model  of  the  heat  balance  between  the 
filament  and  the  surrounding  gas. 

Besides  the  low  frequency  associated  with  the  dynamics 
of  large  scale  buoyancy-driven  vortices,  a  new  set  of  higher 
frequencies  believed  associated  to  the  internal  jet  core  of 
the  flame  has  been  found  in  the  CO-H2  and  N2 -CO-Hs  diffusion 
flames.  Preliminary  results  gave  a  25-30  Hz  range  for  the 
low  frequency  and  300-325  Hz  for  the  high.  However,  a  more 
detailed  Fourier  Transform  frequency  analysis  must  be 
carried  out  to  establish  the  exact  correlation  between  the 
different  harmonics  that  appear  in  the  spectrum. 

The  2-perpendicular  laser-sheets  lighting  technique  has 
been  used  as  a  tool  to  visualize  the  flame  structure 
together  with  the  TFP  technique  to  a  better  understanding  of 


the  combustion  dynamics. 
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Appendix  A: 
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Flame  and  Intensity  Contour  Lines  at  Various 
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Fig.  A-9  Temperature  Profiles  of  N2-CO-H2  Diffusion 
Flame  and  Intensity  Contour  Lines  at  160  mm 
Axial  Location,  experiment  D. 
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Appendix  B:  leffiPfiXSiars  Qi  Ma..-:  gQ  -  Ha  _D  i  fju  3  ign 
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Fig.  B-7  Temperature  Profiles  of  N2-C0-H*  Diffusion 
Flame  and  Intensity  Contour  Lines  at  120  mm 
Axial  Location,  experiment  E. 


90 


9 


Axial  Location,  experiment  E. 


Appendix  C:  Temperature  Prof  ilea  o£  Bif  jus  Aon 

Flams  and  latasaitz  Costoar  Iiinaa  at  Varlsas 
Axial  Locatlona:  Bs  =  2099 
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e  Profiles  of  N2-CO-H8  Diffusion 
Intensity  Contour  Lines  at  0  mm 
Axial  Location,  experiment  A. 
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